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Abstract 

Secondary electron cascades can induce significant ionisation in condensed matter 
due to electron-atom collisions. This is of interest in the context of diffraction and 
imaging using X-rays, where radiation damage is the main limiting factor for achiev- 
ing high resolution data. Here we present new results on electron-induced damage 
on liquid water and ice, from the simulation of Auger electron cascades. We have 
compared our theoretical estimations to the available experimental data on elastic 
and inelastic electron-molecule interactions for water and found the theoretical re- 
sults for elastic cross sections to be in very good agreement with experiment. As a 
result of the cascade we find that the average number of secondary electrons after 
100 fs in ice is about 25, slightly higher than in water, where it is about 20. The 
difference in damage between ice and water is discussed in the context of sample 
handling for biomolecular systems. 
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1 Introduction 



High energy X-rays damage matter, mainly through the photoelectric effect. 
A photon is absorbed by an atom and a photoelectron is emitted from the 
inner (K) shell of the atom. In biological molecules, the remaining hollow ion 
might relax through Auger decay, in which an electron from a higher level 
falls down and a further electron is emitted with an energy between 250 eV 
and 2 keV [1]. The basic atomic physics for these processes is well understood 
and cross sections for photoionization have been tabulated [2,3]. Furthermore, 
Auger line widths have been measured giving K-hole life times, of e.g. 11.1 
fs (C), 9.3 fs (N), 6.6 fs (O) and 1.3 fs (S) [4]. During the photoionization 
process, the photoelectron may interact with electrons in the valence shell. 
For elements of biological significance this may lead to the so called shake-up 
and shake-off effects [5] in which a further low energy electron (10-100 eV) 
is emitted. Semi-empirical quantum calculations were used to estimate that 
such shake-off ionizations occur in 10 - 30% events [6]. 

The advent of X-ray free electron lasers (FEL) will enable a whole range of new 
experiments in physics, chemistry and biology within a few years [7,8,9]. In the 
category of biological applications, we envision performing structural studies 
on large biomolecules, biomolecular aggregates or nanocrystals. Molecular dy- 
namics simulations of a protein molecule in a FEL beam are encouraging as to 
the feasibility of such experiments [10]. Preliminary studies with a simple hy- 
drodynamic model have suggested that [11] a layer of water or ice may delay 
the eventual Coulomb explosion of a biomolecule, such as a protein, in a FEL 
beam [10]. The reason for this is that water molecules at the surface supply 
electrons to neutralize the charging core of the droplet, where the protein is 
expected to be situated, while protons, carrying the positive charge, are ex- 
pelled from the water shell. This may be advantageous in future experiments 
if the sample is held in water droplets. 

In the aftermath of a photoionization event in a liquid or solid, further elec- 
trons will be liberated through collisions between the Auger electron and the 
atoms, leading to a cascade of ionizations. Recent studies [12,13] have made 
a critical assessment of the total amount of electrons due to an Auger electron 
of 250 eV (carbon) 2 . It was predicted that in diamond, between 6 and 13 elec- 
trons will be released, while for amorphous carbon the number is somewhat 
higher due to the absence of a band gap [12]. Here we present results on water, 
both in liquid and solid state. The understanding of Auger electron cascades 
in water is essential for modeling the shielding effects water molecules will 
have on e.g. a protein sample. Furthermore, water is less dense than diamond 



2 Although the actual Auger energy for carbon in diamond is 265 eV [14], we have 
used the value of 250 eV for comparison with earlier work. 
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and therefore more comparable to biomolecules. The Auger electron energy for 
oxygen in a water molecule is 507.9 eV [15], providing twice as much energy 
to the electron than carbon. 



2 Theory 

To calculate the cross sections for scattering of electrons on atoms used in the 
molecular dynamics model, we have used a solid state approach, assuming a 
neutral crystal, or, in the case of liquid water, a neutral cluster with periodic 
boundary conditions. The same approach has been used earlier by Ziaja et al. 
[12,13]. Under this assumption, there are three different interactions that a 
free electron in the medium can undergo: (i) elastic collision with atoms, (ii) 
inelastic collision with atoms or (iii) recombination with atoms or ions. The 
probability of electron-ion recombination is low on the time scale considered 
here (t ~ 100 fs) and therefore neglected [16]. Inelastic collision of electrons 
with ions are not considered, since the number of ionizations caused by a 
photoelectron or single primary Auger electron is much smaller than the total 
number of atoms in the sample. Scattering of electrons on phonons and the 
vibrational cross section are not considered here, though they become relevant 
for low energies (E < 10 eV) [17]. 

Elastic scattering 

The cross section for the elastic scattering of electrons on atoms is calculated 
using programs from the Barbieri/van Hove Phase Shift package [18]. The 
calculations are based on the partial wave expansion technique, using the 
muffin-tin potential approximation. The expression for the elastic cross section 
as a function of the electron energy E, is 



where 5i is the phase shift of each partial wave /. For further details on the 
elastic cross section calculations see reference [13] and references therein. 

The energy of an elastically scattered electron is conserved, and the only ef- 
fect the collision has on the system is that the direction of the electron ve- 
locity changes. For the calculations of the cross sections in liquid water we 
have used a unit cell with hydrogen bonded water molecules at liquid den- 
sity (0.9965 g/cm 3 ), while for ice calculations we have used hexagonal ice, 1^ 
(density 0.9197 g/cm 3 ). The results of our calculations are shown in Fig. 1 
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together with results for diamond [13]. The results compare very well with the 
experimental data for electron scattering on water molecules. 



Inelastic scattering 

A correct treatment of the inelastic scattering is more complicated and, to 
our knowledge, there is no complete method to describe this. We followed 
the approach presented in our previous work [13]. We have used two different 
optical data models, the Ashley model [21,22] and the Tanuma, Powell and 
Penn model (TPP2) [23,24], together with the Lindhard dielectric function 
e(q,u) [25,26,27] for the calculations of the inelastic cross section 



^inelastic 




were hq± = \/2m e (VE ± \JE — huj . The above cross section depends on the 
knowledge of the dielectric function, which is only known for photons, i.e. 
for momentum transfer equal to zero. The two different models, Ashley and 
TPP2, present different ways to empirically compensate for that. The Ash- 
ley approach includes the exchange between the incoming electron and the 
electrons in the crystal, in accordance with the non relativistic M0ller cross 
section [28]. The TPP2 model is developed for calculating the differential in- 
elastic mean free path of electrons in a solid. A more thorough discussion of 
the two models can be found in ref. [13]. From the absorption coefficient for 
liquid water [29], we deduced the dielectric function e(0,u) for zero momen- 
tum transfer. For ice we used experimental data [30] for energies up to 30 eV, 
and for higher energies we have extrapolated the values assuming the shape 
of the curve to be the same as for liquid water (Fig. 2). 

An electron that scatters inelastically will change its path and loose energy to 
the lattice and, if the lost energy is high enough to excite a bound electron from 
the valence band to the conduction band, to the excited electron. For simplicity 
we assumed that the energy lost to the lattice is absorbed into the system and 
in the present study we have not included any treatment for holes. The band 
gaps of water and ice we have used for the calculations were E^ er = 8.7 
eV [31] and E^ p = 7.8 eV [32]. For water the Fermi level E™ ter = 14.85 
eV is taken from experimental data [31], and for ice Ep e = 11.74 eV was 
calculated using the free electron approximation with the density taken from 
reference [33]. The calculations were performed at temperature T lcc = 270 K 
and T water = 300 K. The inelastic cross sections for ice and water are shown 
in Fig. 3 together with those from diamond, taken from reference [13]. 
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From the wealth of experimental data for electron impact ionization of atoms 
and molecules, models and recommended cross sections have emerged for dif- 
ferent targets. We have compared the theoretical estimates for inelastic scat- 
tering of electrons in solids (diamond, ice), with the recommended cross sec- 
tions for carbon [34] and a successful model that describes ionization of water 
molecules [35] (Fig. 4). The former represent a fit to the atomic data for several 
light atoms and ions [34], while the latter, known as the Binary- Encounter- 
Bethe (BEB) model, provides analytical formulae for the ionization cross sec- 
tion of atomic and molecular orbitals by summing up all contributions. The 
BEB model was found to reproduce experimental data for small atoms and a 
variety of molecules very well [35]. 

We find that, in the case of diamond, the inelastic cross section is lower than 
the ionization of single atoms, reflecting the nonlocalized nature of the elec- 
trons in the solid. This difference could prove significant for the estimation of 
electron interaction in solids, if the solid loses its crystalline structure (due to 
Coulomb explosion while exposed to a FEL beam). For ice and water, how- 
ever, we find a smaller difference between the ionization cross sections [35] and 
our theoretical estimates (Fig. 4), and with the opposite sign. 



3 Methods 



Classical simulations of the electron trajectories were performed, where New- 
ton's equations of motion were solved using a Leap- Frog integration scheme [36] 
with a time step of At = 1 attoseconds. The lattice, be it diamond, water or 
ice, was modeled implicitly through the cross sections for collisions. Elastic 
as well as inelastic collisions were modeled as stochastic processes, where the 
probability P{E) for an event was determined from: 

P{E ) = ^a(E) ■ \v e (E)\ ■ At, (3) 



where p is the density of the material, N av is Avogadro's number, m is the 
mass of the atom (C in diamond), respectively water molecule, and v e is the 
velocity of the electron. At was chosen such that P <C 1 in all cases. At each 
time step in the simulation the following processes occur: 
(1) An electron with a certain kinetic energy E^ travels through the sample. 
The scattering probability (both inelastic and elastic) P(E k ) is compared to 
a random number < r < 1. If no scattering occurs, the velocity is kept un- 
changed. 

(2a) If the electron scatters elastically it changes its momentum according to 
a probability function derived from the differential cross section (the formulae 
in [13] were applied for water and ice). No momentum transfer to the lattice 
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has been considered at this point. 

(2b) In the case of inelastic scattering, the electron will lose energy -Ei os t and 
change its path. If £i ost > E gap an additional electron is liberated from the va- 
lence band with E k = Ei ost — E gap . If E\ ost < E glLp , the energy is transfered to 
the lattice. In both cases the momentum is changed based on the differential 
cross section calculated for water and ice using formulae from [13]. 
(3) The electron position is integrated. 

The spatial electron dynamics program is part of the GROMACS software 
package [37]. 



4 Results and Discussion 

The classical dynamics simulations of the electrons produced the same results 
for diamond as our earlier Monte Carlo simulations [12], showing that the im- 
plementation of the algorithm is correct, and that the integration time step 
was small enough. In principle, the model could be enhanced by adding electro- 
static interactions between electrons and ions [38]. It is however not trivial to 
model the interactions between electrons and ions classically without severe 
approximations. Since we have only simulated single Auger electron trajec- 
tories, and the cross sections were determined for neutral systems, we have 
deemed it unnecessary to add the Coulomb term. Finally, we note that correct 
treatment of electron/ion or electron/hole systems can be performed using e.g. 
time-dependent density functional calculations, but it is complicated, and not 
computationally feasible for macroscopic systems [39]. 



Ionization rates 

We find that an Auger electron of 500 eV will produce typically between 5 and 
10 secondary electrons in the first femtosecond in ice and water, a result which 
is comparable with the average number of ionization events in diamond from 
a 250 eV Auger electron. In ice and water, however, saturation appears much 
later, after about 100 fs, when the average number of ionizations is around 20 
for liquid water and 25 for ice. These results are shown in Fig. 5, which not 
only compares the ionization rates for ice, water and diamond, but also shows 
that the simulation of diamond using classical dynamics reproduces previous 
results obtained with a Monte Carlo method well [12]. We note that, although 
the cross sections for ice and water are rather similar, the difference in the 
number of ionizations in the Auger electron cascade is determined mainly by 
the the approximately 1 eV difference in the band gap between liquid water 
and ice. 
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Energy distribution 



The total energy of the electrons shows a substantial drop during the first 
10 femtoseconds, with the rest of the energy being transfered to the lattice 
(Fig. 6). The electron gas is cooling down very rapidly, due to the large number 
of ionizations, and the average temperature of the electrons kT e \ is similar for 
ice and water. The energy lost to the lattice is more than 50% of the initial 
Auger electron energy. For comparison, the results obtained for diamond are 
also shown, consistent with the previous calculations [12]. 

The electron energy distribution as a function of time in the cascade (Fig. 7) is 
particularly interesting. In the early stages of the cascade (of the order of 0.1 
fs), the electron population is divided in two sets, initial Auger electrons which 
are highly energetic (450-500 eV) and the first secondary electrons of lower en- 
ergy. At the critical time of about 1 fs, the population of high energy electrons 
fades away, leaving a spread of electrons with intermediate and low energies. 
Subsequent collisions will cool down the electron gas very fast, and after 10 fs 
all the electrons will have energies less than 20 eV. At such low energies (be- 
low 10 eV) and long time scales (100-1000 fs), new processes become relevant, 
such as electron-ion recombination or electron-phonon scattering, which are 
not yet included in our model. 

Spatial distribution 

We have performed an analysis of spatial distribution of electrons in the cas- 
cade in our classical dynamics model for ice, water and diamond. Fig. 8 shows 
the average gyration radius of the electron cloud for ice, water and diamond. 
The spatial extension of the cloud in ice and water is three times larger than 
the one in diamond, due to the considerably lower density and the correspond- 
ing longer mean free path an electrons in water/ice. 



5 Conclusion 

An understanding of the behaviour of water and ice exposed to an intense X- 
ray beam is important for planned biology experiments with FELs. Biological 
macromolecules carry about 30-60% water with them as structural water and 
water bound as a surface layer. In the light of results from a hydrodynamic 
model [11], it was suggested that embedding a sample in water or ice may 
be beneficial in the sense that the eventual Coulomb explosion [10] could be 
delayed by the presence of the water or ice layer. Our results show that in 
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hexagonal ice, Ih, a slightly higher number of secondary electrons are liber- 
ated as compared to water (Fig. 5), despite the slightly lower temperature. 
Furthermore, the electrons are spread over a larger volume (Fig. 8). Although 
the differences are small, it seems therefore that water is the preferred medium 
for encapsulating biomolecules when studied in FEL beams. If the sample is 
embedded in vitrious ice, which could be true in case of a rapid freezing, we 
expect the behaviour to be unaffected compared to water. 
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Figure captions 

(1) Elastic cross sections for scattering of electrons on ice, water and dia- 
mond. Lines represent the scattering cross section per scattering center, 
derived using the Barbieri/van Hove Phase Shift method [18], for ice (full 
line) and water (dashed line). For comparison, the cross section for di- 
amond taken from [13] is shown (dotted line). The experimental points 
represent the elastic scattering of electrons on water molecules as mea- 
sured in [19,20]. 

(2) The imaginary part of the inverse dielectric function, Im(— l/e(0, E)), for 
water and ice plotted as a function of the incoming photon energy. 

(3) Inelastic cross section for electron scattering on ice and water as a function 
of the electron energy. Lines represent the cross section per scattering 
center, calculated using two optical data models, the Ashley model [21,22] 
and the Tanuma, Powell and Penn model (TPP2) [23,24]. For comparison, 
the cross sections for diamond taken from [13] are also shown. 

(4) Comparison between the inelastic cross section for electron scattering 
on solids (ice and diamond) and the ionization cross section for elec- 
tron scattering on single atoms and molecules. Lines represent the cal- 
culated cross section using TPP2 model [23,24] for diamond (solid line) 
and ice (dashed line). Points represent the recommended cross section for 
carbon atoms [34] (solid points), and the calculated cross section Binary- 
Encounter-Bethe (BEB) model [35] for ionization of water molecules (cir- 
cles). 

(5) Average number of secondary electrons as a function of time, from an 
initial Auger electron (at time t=0 fs) in ice (solid line), water (dashed 
line) and diamond (dotted line). Lines correspond to calculations using 
Ashley's model [21,22]. 

(6) Average energy of the free electrons in an electron cascade coming from 
an initial Auger electron in ice and water as a function of time (thin 
lines). The rest of the Auger electron energy is transfered to the lattice. 
Average temperature kT e \ of the electron gas versus time (thick lines) for 
water and diamond. 

(7) Histogram of the energy distribution N(E,t)/N(t) among electrons at 
different time steps in ice. At a fixed time, the histogram is normalized 
to the number of electrons in the cascade N(t). Results are shown for an 
average cascade produced by an Auger electron in ice. 

(8) The gyration radius (average distance of electrons from the cloud center 
of mass) of the electron cascade as a function of time, for ice, water and 
diamond. 
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